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Two-dimensional Packing in Prolate Granular Materials 

K. Stokely, A. Diacou, and Scott V. Franklir|3 

Dept. of Physics, Rochester Institute of Technology 
(Dated: February 1, 2008) 

We investigate the two-dimensional packing of extremely prolate (aspect ratio a = L/D > 10) 
granular materials, comparing experiments with Monte-Carlo simulations. In experimental piles of 
particles with aspect ratio a = 12 we find the average packing fraction to be 0.68±0.03. Both exper- 
imental and simulated piles contain a large number of horizontal particles, and particle alignment 
is quantified by an orientational order correlation function. In both simulation and experiment the 
correlation between particle orientation decays after a distance of two particle lengths. It is possible 
to identify voids in the pile with sizes ranging over two orders of magnitude. The experimental 
void distribution function is a power law with exponent — /3 = —2.37 ± 0.05. Void distributions in 
simulated piles do not decay as a power law, but do show a broad tail. We extend the simulation to 
investigate the scaling at very large aspect ratios. A geometric argument predicts the pile number 
density to scale as a -2 . Simulations do indeed scale this way, but particle alignment complicates 
the picture, and the actual number densities are quite a bit larger than predicted. 
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I. INTRODUCTION 



One of the more striking features of piles of very pro- 
late granular materials (large aspect ratio a = L/D) is 
the connected network that forms at comparatively low 
packing fractions. The formation of this network is of- 
ten commercially undesirable. LCD screens, for example, 
cannot function if the molecules are entangled, and lum- 
ber floating down a river stops when the logs jam. There 
are, however, practical applications for such "jammed" 
networks. Piles of large aspect-ratio materials are ex- 
tremely rigid, even at low packing fractions, and have a 
high strength:weight ratio. At the extremely small scale, 
networks of carbon nanotubes are a possible mechanism 
for conducting energy to and from nano-devices pj. 

Little is known about even basic characteristics of piles 
formed from rod-like particles, most research on non- 
spherical particles, whether in two [J 0, IE IE S 13 or 
three |2( dimensions, being limited to a < 5. The rigidity 
of such piles is due to particle entanglement, with particle 
rotation extremely constrained. The statistics of particle 
orientations which determine these constraints, however, 
is not known. While it seems obvious that particles will 
align, in fact two-dimensional piles contain a number of 
orthogonal particles that create large voids which domi- 
nate the pile landscape (see Fig.^l. The only work above 
a ~ 10 we are aware of is that of Philipse [E Eli wn0 
formed three-dimensional piles of copper wire of aspect 
ratios ranging from 5 to 77 and explained the 1 /a scaling 
of the volume fraction with a simple geometric model. As 
the particles' aspect ratio increased, they could no longer 
be poured from their initial container, and tended to fall 
out as a solid "plug" . The cause of this transition, which 
occurs at a ~ 35, is not known. 




FIG. 1: Top: 2-d pile of a = 12 acrylic rods backlit with 
fluorescent lights. The middle third of each particle appears 
bright. Bottom: Simulated pile of a = 12 particles. Both 
piles show particles aligning and a wide distribution of void 
sizes. 



II. EXPERIMENT AND SIMULATION 

To form prolate particles, acrylic rods (diameter D = 
0.16 cm) were cut to a length L = 1.9 cm (a w 12) and 
constrained between two Plexiglas sheets separated by a 
spacer 1.25 particle diameters thick. The uniform spac- 
ing throughout the plates prevents particles from overlap- 
ping; piles are effectively 2-dimensional. Thicker spacers 
result in overlapping particles which are pinched between 
the plates and immobile. Whereas Philipse observed a 
qualitative increase in pile rigidity in three dimensions 
to occur at about an aspect ratio of 35, we believe that 
in two dimensions this occurs at much lower aspect ratios. 
We have observed 2-d piles of particles with aspect ratio 
10 to have stable angles of repose of 90° or greater (see 
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FIG. 2: 2-dimensional pile of a = 10 particles showing an 
angle of repose of 90° . We believe that the ability to support 
angles of 90° or greater is analogous with Philipse's observa- 
tion of solid plugs in three dimensions. 



solid plug observed in 3-d. 

Piles arc formed by distributing particles at random 
on one Plexiglas sheet, attaching the second sheet, and 
slowly rotating the system to vertical. The initial dis- 
tribution is not truly random; local orientational corre- 
lations exist as neighboring particles are forced to be 
aligned (or else they would overlap). This could be 
avoided in principle by reducing the number of parti- 
cles on the plate; in practice, however, this would be 
prohibitively slow. Additional particles are prepared in a 
similar manner in an identical cell which is used as a fun- 
nel to pour particles onto the pile. Piles formed this way 
are 53 cm wide, typically 25 cm high and contain about 
2000 particles. A picture of a pile is shown in Fig. 
(top). The piles are backlit with fluorescent lights. The 
cylindrical rods act as lenses, displaying a thin bright line 
throughout the middle of each particle. We have written 
software that identifies connected, collinear bright pixels 
in a picture and extracts the particle location and orien- 
tation. Data reported involve averaging over 19 separate 
piles. Despite the less-than-ideal preparation, piles were 
statistically consistent; packing fractions varied by ~ 5% 
and void distribution and orientational order functions 
were similarly reproducible. 

Buchalter and Bradley^, 0] developed a Monte-Carlo 
simulation for ellipsoidal particles. We adapted this for 
cylindrical particles and extended the aspect ratio by two 
orders of magnitude (a max = 1000). Particles move along 
the nodes of a discrete lattice (N x N, with N ~ 10L) 
and can rotate freely. Particles are initially placed at 
random locations on the lattice and given random orien- 
tations (never being allowed to overlap with other parti- 
cles). A single particle is chosen at random and moved 
along a randomly generated displacement /rotation path. 
The only constraint on the motion is that particles can- 
not move upwards or overlap with other particles. The 
maximum possible distance a particle can move in one 
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particle is placed at its last allowable position and a new 
particle is then chosen for an attempted move. At any 
given time, only one particle is in motion. The process 
repeats until the potential energy (the sum of the parti- 
cle heights) remains constant for 5000 time steps, each 
particle unable to move for, on average, 10 attempted 
moves. A new group of particles is then placed above the 
formed pile and allowed to settle. All piles are at least 7 
particle lengths high and we have checked to ensure that 
additional pourings do not appreciably change the pile's 
statistics. A sample pile is shown in Fig. ^ (bottom). 
The length of the particles, constant through any one 
pile, is varied from 10 to 1000. Results for a given aspect 
ratio are averaged over five piles; additional piles do not 
change the statistics. Additional details about the sim- 
ulation's validity, including a discussion comparing the 
simulation in the limit as the aspect ratio goes to 1 with 
experimental findings, can be found in 

III. RESULTS 
A. Global Pile Characteristics 

The range of packing fractions achievable with 2-d 
disks under gravitational forces is quite narrow. The 
upper and lower limits are given by hexagonal (<j>hcp = 
tt/(2v / 3) ~ 0.907) and orthgonal (<f> ocp = tt/4 « 0.785) 
close packing respectively, with a random close packing 
vale of (f> rcp w 0.82 llj. We find the average packing 
fraction of rods with a = 12 to be = 0.68 ± 0.03. The 
lowest measured value was 0.63; the largest 0.72. 

The orientational order parameter Q = (cos(2#i)) is 
used to characterize the angular distribution of particles. 
6i is the angle with respect to the horizontal of the i-th 
particle and the average is over all particles. Q takes 
values ranging from —1 (all vertical) to +1 (all horizon- 
tal), with Q = indicating an isotropic distribution of 
angles (or all angles equal to 45°). For the experimental 
piles Q — 0.33 ± 0.05; simulated piles have similar values 
regardless of aspect ratio. 

Figure|3|shows the distribution of particle angles in ex- 
perimental and simulated piles. Both have a peak around 
8 = 0, indicating a preference for horizontal orientation; 
this preference is stronger in the simulated pile. To deter- 
mine whether this ordering was caused by the flat bottom 
boundary, we calculated the orientational order parame- 
ter for all particles whose centers of mass are at height 
h. We denote this height dependent value as Qh{h). A 
plot of Qh(h) vs. height for experimental piles is shown 
in Fig. 2| Qh(h = 0) is 1, as particles on the floor must 
be horizontal. At a height of h = L, however, Qh{h) has 
already decayed significantly. For heights greater than a 
particle length, Qh(h) fluctuates about an average value. 
From this we infer that the bottom boundary's influence 
on particle orientation does not extend beyond one parti- 
cle length. Simulations show a similar asymptotic value 

for C)iA h) anrl it should be noted that the bnttnm bmmrl- 
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FIG. 5: Void distribution function V(A) as a function of void 
FIG. 3: Angular distribution of experimental (•) and simu- area scaled by particle length L squared. V(A) is expressed 
lated (*) particles. Both show a peak about horizontal align- as a percentage of the total number of voids. Voids in both 
ment (9 — 0), the experiment having the broader distribution. experiment (•) and simulation (*, with a = 20) decay as a 

A' 13 with p = 2.37 ± 0.05. This means that smaller voids 
occupy a greater total area than the larger, rarer, voids. 



V 



0.8 



0.6 



Q h (h) 



0.4 



0.2 



* ****** 



* * * 



1 2 

Height / Particle Length 

FIG. 4: Orientational order parameter Qh(h) as function of 
height in pile. The reaching of an asymptotic value so quickly 
(for h ~ L) indicates that the bottom boundary is not sig- 
nificantly influencing the orientation of particles high up in 
the pile. Simulations, which do not have a horizontal bottom 
boundary but rather allow the bottom particles to assume 
any angle, show a similar average Qh(h) at all heights. 



B. Distribution of Voids 

The appearance of both simulated and experimental 
piles are dominated by large, but rare, voids. From the 
images we find the number of voids as a function of void 
area A. The void distribution functions V(A) from ex- 
perimental (•) and simulated (*) piles are plotted vs. 
void area A in Fig. [5] As Fig.[S]shows, experimental void 
sizes vary by over two decades. The experimental data 
are well-fit by a power law A~ 2 - 37±005 (straight line in 
Fig. 0. Simulated piles show a similar decay, although 
the function does not seem to follow a power law. We did 
not notice any significant dependence of the void distri- 
bution function on the particle aspect ratio in simulated 
piles, although this warrants further study. 

There are several intriguing consequences to the fact 
that the exponent in V(A) is between -2 and -3. First, 
the total area taken up by all voids with size A is 
AV(A) oc A' 137 . As lim A ^ 00 [At/(A)] -> 0, the cumula- 
tive effect of the smaller voids to the pile's area is actually 
larger than that of the larger voids. We also note that 
the total area occupied by all voids J AV(A)dA remains 
finite as, in fact, it must for realistic piles. The total area 
occupied by voids is related to the packing fraction <f> by 
the relation 



ary in the simulations does not impose a horizontal angle 
on the bottom particles. Therefore we believe the ten- 
dency for particles to be horizontal, more pronounced the 
simulation than in the experiment, is a result of gravity 
(in the simulation imposed by the restriction that parti- 
cles cannot move upwards) rather than a boundary condi- 
tion. Experimental piles also have more vertical particles 
than the simulation, a consequence, we believe, of fric- 
tion between particles which is not incorporated in the 

simulation 



AV{A)dA = / A dA = A tot (l - 0) 



where A tot is the total pile area. This relation can be 
used to explore the lower limits for the void size, with- 
out which the integral on the left diverges. This is the 
subject of current research. We also note that the ab- 
sence of an upper limit for void size results in the di- 
vergence of the integral for the mean square void area 
1 4 2 \ = f V(A\A 2 dA 
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C. Neighboring Particle Alignment 

Both experiment and simulation find neighboring par- 
ticles aligning. This is quantified with an orientational 
correlation function 

Q(r) = (cos(2A%)). 

AOij is the difference in angle between particles i and 
j and the average is over all particles whose centers-of- 
mass separation is between r and r + Sr. This func- 
tion, related to an earlier orientational order parame- 
ter |3j, takes values ranging from 1 if particles are 
parallel to -1 if particles are perpendicular. Two parti- 
cles whose centers-of-mass are quite close must be aligned 
and so Q(r — > 0) — > 1. Once particle centers-of-mass 
are separated by more than one particle length L they 
can in principle assume any relative orientation and so 
Q(r — > oo) — > 0. For comparison, we calculate analyt- 
ically the correlation function resulting from a distribu- 
tion where particless assume all allowable angles with 
equal probability. This is the simplest possible order- 
ing, the only constraint being that particles cannot over- 
lap, and is used in simple geometric models for predict- 
ing number density. The allowable angles a particle 
can take with respect to a fixed particle assumed to lie 
along the x axis are found as a function of center-of- 
mass separation r and angle <j) that the line connecting 
the centers-of-mass makes with the x axis. If the mini- 
mum/maximum allowable angles are given by 8 m i n and 
9 max then Q(r) is 

1-2-k rO max (r,<f>) 

Q{r) = d<t> cos (20)d0. 

Figure El shows the Q(r) distribution resulting from the 
analytic (line), experimental (•), and simulated (*) piles. 
Both the experimental and simulated piles show greater 
correlation between neighboring particles, seen in the di- 
vergence from the analytic line for r/L > 0.5 (between 
the dashed lines), and reach an asymptotic value once 
particles are separated by more than two particle lengths. 

The long-range correlation between particles shown in 
Fig. does not represent a long-range influence of one 
particle on another, but rather results from the overall 
preference for particles to be horizontal. This is con- 
firmed by calculating 

Q(r -► oo) = J P(0)P(<j>) cos [2(0 - <f>)]d0d(f>, 

which assumes the particle angles are drawn at ran- 
dom from the distribution P{9) shown in Fig. [21 The 
differences in the simulation and experimental distribu- 
tion functions result in Q c ^ = 0.16 and Q™ 11 = 0.53, 
agreeing quite well with the asymptotic values in Fig. [S] 
Q(r/L) reaches its asymptotic value for both simulation 
and experiment within two particle lengths. This cor- 
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FIG. 6: Orientational correlation function Q(r/L) as a func- 
tion of center-of-mass separation scaled by particle length. 
The solid line represents the correlation resulting from a dis- 
tribution where particles assume all allowable angles with 
equal probability. Both experiment and simulation show en- 
hanced alignment for r/L between 0.5 and 1 (between dashed 
lines) but noticeably different asymptotic values due to the 
different angular distributions from Fig. |3] 




FIG. 7: Orientational correlation function Q(r/L) normalized 
by the asymptotic value Qoa as a function of center-of-mass 
separation scaled by particle length. Simulation results are 
shown for particle lenghts varying over three orders of mag- 
nitude; the correlation length at which the asymptotic value 
is reached is the same for simulations and experiment. 



lengths differing by two orders of magnitude, as shown 
in Fig. [7| The curves in Figure [7| have all been normal- 
ized by their asymptotic value; that is, what is plotted 
is (Q — Qoo)/(l — Qoo), which takes values ranging from 
1 to 0. When thus normalized, all simulated curves lie 
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FIG. 8: Top: Contact number in simulated piles as a func- 
tion of particle aspect ratio, (c) is independent of aspect ratio 
for long particles, and reaches a suitable asymptote, 4.5, com- 
parable withe contact number of disks, as a — > 1. Bottom: 
Orientational order parameter Q as a function of aspect ratio 
for simulation (*) and experiment (•). 



D. Simulation at Large Aspect Ratios 

We now extend the simulation to larger aspect ratios 
and investigate the scaling of various quantities. The 
rigidity of a pile depends on particles in contact, hence 
we calculate the contact number (c). Fig. |Hftop) shows 
that (c) reaches an asymptotic value of ~ 3.2 by aspect 
ratio 50. This may seem counter-intuitive, as longer par- 
ticles in principle can be in contact with more neighbors. 
Orientations that maximize (c), however, are quite rare 
and the length-independence of contact number is due to 
the tendency of neighboring particles to align and screen 
one another from other particles. As particle aspect ra- 
tio decreases, and the particles become more circular, 
this screening effect diminishes and the contact number 
increases. Packings of circular disks, for example, show 
a contact number between 4 (orthogonal close packing) 
and 6 (hexagonal close packing). 

Recall that the pile's orientational order is character- 
ized by the order parameter Q = (cos(2f?;)) where Oi is 
the angle the ith. particle makes with the horizontal and 
Q is averaged over all particles. As shown in Fig. 03 Q of 
simulated piles decreases as the particle length increases 
reaching an asymptote of 0.3. Q for experimental piles 
of particles with aspect ratio a — 12 is 0.33 ± 0.05, com- 
parable with that of of simulations. 



E. Phenomenology 

Philipse gave a simple geometric argument, the Ran- 
dom Contact Model 9] (RCM)to explain the low pack- 
ing fractions of three-dimensional piles. We now apply 
his logic to our two-dimensional piles and show the dis- 
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FIG. 9: Number density as a function of particle length. 
As predicted by the simple geometric model described in the 
text, the number density appears to scale as L~ 2 for large L. 
The actual densities however, are larger than that predicted 
by the model, a consequence of particle alignment. 



scribed above. 

The existence of one particle excludes a fraction of the 
possible orientations, called the excluded area A exc \, that 
can be assumed by a second particle. If we assume a 
connected network, where all particles are in contact with 
(on average) (c) neighbors and that particles assume all 
allowable orientations with equal probability, then the 
average number density will be 



(N) 



Acxc 



The factor of 2 accounts for the fact that each contact 
involves two particles. We have already shown, however, 
that the assumption that contacts are uncorrelated, is 
not satisfied. Balberg^ has calculated the excluded 
area of a stick with length L and width W; to first order 
A:xci = (2/7r)£ 2 - With (c) = 3.2, the number density as 
a function of aspect ratio a is predicted to be 



N(L) = 



2(c) 



(2L 2 /ir) 



= CL~ 



with C = 10. Fig. ©shows that N(L) does indeed fall off 
as L~ 2 for large L. The constant, however, is larger than 
that predicted by the RCM (flat line in Fig. EJ). Piles 
are therefore more dense than predicted, implying that 
the excluded area of particles is about 33% less than that 
in an isotropic distribution. This is a result of particle 
alignment, seen earlier in Fig. [H] We also note that the 
scaling as or 2 is realized only for the largest of aspect ra- 
tios, while the constancy of contact number occurs much 
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IV. CONCLUSIONS 

We have presented the first quantitative characteriza- 
tion of two-dimensional piles formed from prolate (a > 
10) granular materials finding, for example, the pack- 
ing fraction of particles with aspect ratio a — 12 to be 
0.68 ± 0.03. Particles separated by less than two parti- 
cles lengths show a greater orientational correlation than 
would be found in a random pile; particles separated by 
more than two lengths are uncorrelated except for the 
general preference for horizontal alignment imposed by 
gravity. The void distribution function in experimental 



piles obeys a power law with exponent —/3 — 2.37 ± 0.05; 
Monte-Carlo simulations show similar angular correla- 
tions and void distribution functions. Simulations have a 
greater number of horizontal particles, however, and thus 
produce piles with larger number densities than found in 
both experiment and simple geometric models. 

We thank E. F. Rcdish for first questioning the char- 
acteristics of pickup sticks and John C. Crocker for call- 
ing attention to the work of Philipse. Eric R. Weeks and 
L. S. Mcichlc have provided invaluable advice throughout 
this project. Saul Lapidus and Peter Gee were involved 
in much of the original setup of the experiment. 
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